Adiabatic stabilization of excitons in intense terahertz laser 
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High-precise calculation of near-infrared absorption and transient spectra in bulk semiconductors 
irradiated by an intense terahertz (THz) laser shows that, the ionization rate of the ground-state 
exciton, probed through the dynamic Fano resonance, may decrease with the increase of the THz 
laser intensity. This counterintuitive effect indicates the excitons are stabilized against the field- 
ionization. The much lower Rydberg energy and "atomic unit" of laser intensity for excitons, and the 
possibility of creating excitons from the "vacuum state" allow observing the excitonic stabilization 
in experiments, in contrast to the case of atomic stabilization. 



PACS numbers: 71.35.Cc, 78.20.Jq, 32.80.Fb, 42.50.Md 



In accordance to Einstein's theory of the photoelectric 
effect, an electron can be stripped from its atom by light 
beams with sufficiently high frequency. As the intensity 
of a light beam is increased, the stripping probability 
of the electron increases owing to the increasing number 
of photon impacts. Surprisingly, since the late 1980's, 
theoretical investigations and numerical simulations 
have predicted that when the electric field, associated 
with a laser of sufficiently high intensity and frequency, 
approaches or exceeds the electrostatic field between the 
electron and the ion in an atom, the wave function of 
the irradiated atom may be distorted adiabatically into 
a distribution with two well separated peaks. The peak 
spacing increases with increasing the field intensity, thus 
the atomic electron spends more time far away from the 
nucleus, and the ionization rate slows down dramatically 
till almost totally suppressed. In other words, the adia- 
batic stabilization of the atom occurs. 

The observation of this atomic stabilization effect is, 
however, extremely difficult. Firstly, for the atomic 
ground state, the high-frequency and high-intensity con- 
dition requires an extreme ultraviolet laser with intensity 
larger than an atomic unit (~ 3.51 x 10 16 W/cm 2 ), which 
seems not available in the near future. Secondly, and 
more vitally, sufficiently slow turn-on is demanded for 
the laser pulse to adiabatically drive the atomic ground 
state into a stable dressed state; during the rise time of 
the super-intense laser pulse, however, substantial ion- 
ization has already been inevitable |Q, which prevents 
the effect from observation. Some researchers even ar- 
gued that it was impossible to observe such an adiabatic 
stabilization of atoms ^ n ^ ac ^' except for a few dis- 

putable indications for the stabilization of high Rydberg 
states Jl(i[| , so far there is no unambiguous experimental 
evidence for the atomic stabilization. 

It is compelling to settle the debate whether this effect 
exists or not. Apart from the pure curiosity in basic re- 
search, the atomic stabilization is also of importance in 
potential application, such as the high harmonic genera- 



tion [ pd| and quantum computation |l2| ] , as it can quench 
the undesired ionization and dephasing caused by the in- 
tense laser that is used to drive electrons or to manipulate 
a qubit in an atom. 

In this Letter, we report on the adiabatic stabiliza- 
tion of hydrogen-atom- like excitons in semiconductors. 
The excitonic stabilization (XS) , in contrast to its atomic 
counterpart, is definitely an observable effect, owing to 
the basic characteristic of excitons that they are created 
from "vacuum" of carriers by the interband optical exci- 
tation. In our study, the dressed exciton states are pre- 
pared by weak near-infrared (NIR) laser pulses in semi- 
conductors in the presence of a quasi-continuous tera- 
hertz (THz) field, as in most recent experiments on the 
dynamical Franz-Keldysh effect |]l3|-|l5||. The ionization 
rate or lifetime of the dressed excitons can be directly ex- 
tracted by monitoring either the linewidth of absorption 
peaks or the decay time of transient signals. Thus the 
"turn-on" problem is safely circumvented as long as the 
pulsed interband excitation and the following dephasing 
process of excitons are well covered by the flattop region 
of the THz laser pulse. 

Moreover, the present laboratory gear is already ready 
for observing the XS effect. Due to the small effec- 
tive mass \x and large dielectric constant e, an exciton 
has much smaller binding energy Eb than the hydrogen 
atom does, and the "excitonic unit" for laser intensity 
(cx /i 4 /e 5 5 ) in a typical semiconductor like GaAs is about 
10 orders of magnitude smaller than the atomic unit. So, 
the "high-intensity and high-frequency" requirement for 
the XS can readily be fulfilled by the free-electron lasers 
operating with MW/cm 2 power and THz frequency [fl6l| . 

Our investigation is based on the high-precise calcula- 
tion of interband optical spectra, which includes properly 
the Coulomb interaction, the nonperturbative ac-field, 
and the contribution of continuum states. To avoid un- 
necessary complexity, the excitonic state is treated within 
the simple effective-mass approximation. Then the mo- 
tion equation for the relative motion of an electron-hole 
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pair is 



iWt^T, t) = (H- i 72 )^(r, t) + x (t)e- i;not 5(r), (1) 



H = E„ — --V 2 

9 2^ r iirer 



+ eFzcos(ujt), (2) 



where E g is the band gap, F is the field strength of the 
THz laser associated with frequency uj and linearly polar- 
ized in the z direction, x(t) denotes the NIR pulse exci- 
tation centered at Qq, and 72 is an interband dephasing- 
ratc due to phonon scattering. The THz laser is assumed 
in a continuous wave (cw) form, because its microsec- 
ond duration Jl6| is much longer than the picosecond 
interband excitation and dephasing process. This model 
system is mathematically equivalent to a hydrogen atom 
in an intense laser except that the electron-hole pair is 
generated by the NIR laser. 

Without the Coulomb coupling, the eigensolution of 
the time-dependent Hamiltonian ( Eq. (||)) is the Volkov 
state u%, which, as the basis set, reads 
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where m and k, associated with quasienergy £kr, 
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Up + mhuj, denote the sideband index and the 

is the effective 



wavevector, respectively. Here U p 
ponderomotive potential, dp 



-^2 is the classical ex- 
cursion amplitude of an electron in the ac-electric field, 
and |k) is the Bloch-state with accelerating quasimomen- 
tum k — |^zsin(o;i). With the excitonic Rydberg unit 
adopted, the matrix elements of the Coulomb potential 
read 



qm+n,km 



-2 — 2 -n t / \ 
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where J n (x) denotes the nth-order Bessel function of the 
first kind. The intra-sideband (n = 0) interaction is di- 
rectly related to the Kramers-Hennerberg (KH) poten- 
tial, the time-average of the Coulomb potential in the 
coordinate frame rested on the quivering electron |jl8| . 
The KH potential has logarithmic singularity in the seg- 
ment between the turning points ±za^ and r" 1 / 2 singu- 
larity at the ends. Consequently, the ground state |3>) is 
stretched along the segment of singularity and becomes 
dichotomous for super-intense field (ap 3> ag, the exci- 
ton Bohr radius) |l[ . 

The spectrum of the system exhibits the sideband 
structure. As illustrated in the lower inset in Fig. |], each 
sideband consists of several discrete states and a contin- 
uum, and the discrete states are embedded in the ioniza- 
tion continua of lower sidebands. As predicted very re- 
cently in an ac-driven biased semiconductor superlattice 
p9[ , the quantum interference between a discrete state 
of quasienergy excitons and the degenerate continuum 
of a sideband can lead to the dynamic Fano resonance 
(DFR), which manifests itself as broadened asymmetric 



lineshape in absorption spectra and as intrinsic decay in 
transient four-wave mixing signals. Thus, the ionization 
rate of the quasi-bound Floquet-state excitons and the 
XS effect can be studied via the DFR induced by the 
inter-sideband coupling (n ^ 0) in the present system. 

To calculate the time-dependent wavefunction, we nu- 
merically integrate Eq. (|]J) in the cylindrical polar coor- 
dinates by extending the space-time difference method 
p0| to a time-periodic system, with the initial condi- 
tion as ip(r, —00) = 0. The out-going waves are ab- 
sorbed at the boundary with the mask technique p|. 
The NIR absorption spectrum a(fl) oc 3[P(f2)/x(^) , 
in which P(t) = ip(0, t) gives directly the transient in- 
terband polarization (i.e. coherence). In the calcula- 
tion, huo = 10 meV (~ 2.4 THz), 72 = 1 meV, and the 
NIR pulse is assumed to be of a Gaussian shape, i.e. 
\{t) = Xoexp (— i 2 /2r 2 ). For cw absorption, r is taken 
as 8/w, which is large enough to eliminate the sideband 
overlap. The material parameters take bulk GaAs as an 
example (Eb = 4 meV and as = 10 nm). 

Figure [l] plots the NIR cw-absorption spectra of the 
bulk GaAs driven by the THz-field with various strength. 
The blue-shift of the peak associated with the Is exciton 
state results mainly from the effective ponderomotive en- 
ergy Up and additionally from the decrease of the bind- 
ing energy Eb- The sideband structure of the Floquet- 
state exciton is visible in the spectra, and particularly, 
for F > 24 kV/cm, or the intensity / > 2.3 MW/cm 2 
pjj , the — Ilo sideband can be stronger than the Olu one 
that evolves from the original Is state. There is also 
significant absorption in the band gap. Indeed, the dy- 
namical Franz-Keldysh effect [ p"5[ is properly reproduced 
in this high-precise calculation, and the result can be well 
understood in the Floquet-state picture. 

The exciton ionization rate is probed through the 
broadening and distortion of absorption peaks in the 
NIR-spectra. As displayed in the inset, the asymmetric 
resonance peak is well fitted with the Fano lineshape p^ j 
characterized by the shape parameter q and broadening 
constant 7 

(qj + Q-E g +£ B ) 2 



a{n) = ao + ai ^+(n-E g + £ B r> 

indicating the occurrence of the DFR. By fitting the 
Out peak, the ionization rate T can be extracted via 
7 ~ 72 + r/2. As shown in Fig. |^, 7 (hence the ioniza- 
tion rate) increases with the laser intensity until F « 16 
kV/cm (correspondingly, / w 1.0 MW/cm 2 [^TJ and 
a,F ~ 1.28as). Beyond that critical field strength, the 
ionization rate commences to decrease, and approaches 
zero at F « 28 kV/cm, clearly demonstrating the XS 
effect. 

r can also be obtained analytically in the high- 
frequency approximation [p3| . When the wavefunction 
dichotomy is negligible for medium of (say clf < 3as), 
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a compact result for T of the ground state 
derived as 



231 can be 
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J^(aFk n x)dx, 



where fc 2 , = nuj — £b- When field is enhanced so that 
approaches or exceeds 1, due to rapid oscillation 
of the Bessel function, the integration over x diminishes. 
Physically, this originates from the destructive interfer- 
ence of the out-going waves scattered from different part 
of the initial state. On the other hand, the stretching 
of the state along the field will reduce the wavefunction 
peak, contributing additionally to the XS. 

Fano interference is one of fundamental channels for 
irreversible decay Q,^5|, so the XS can also be directly 
observed in the time domain by monitoring the transient 
optical signal after a pulsed excitation. The time-resolved 
interband polarization is calculated for several THz-field 
strengths and plotted in Fig. [| The duration of the NIR 
pulse r is then set to be 2/w, and the central frequency 
is chosen for each field strength to be resonant with the 
brightest transition, namely Qq — E g = —4, —2, 2, 5, and 
6 meV for F = 0, 8, 16, 20, and 30 kV/cm, respectively. 
To single out the dephasing induced by the DFR only, 
the reformulated signal V(t) = P(t) exp(j2t) is plotted. 
At finite field, the optical signal roughly oscillates with 
the period of T = 2n/ui, because the states in the KH po- 
tential oscillate periodically in the laboratory frame. The 
modulation beat superimposed on the fast oscillation re- 
sults from the quantum interference between the Floquet 
states with different quasienergy. The dynamic Fano in- 
terference between the discrete state and the continuum 
opens a new channel for interband dephasing, leading to 
the decay of the transient signal. The XS effect is ver- 
ified by the fact that the dephasing rate calculated for 
F > 16 kV/cm begins to decrease, consistent with the 
variation trend of the ionization rate as shown in Fig. |^. 
For F = 30 kV/cm (I w 3.6 MW/cm 2 Q), the decay 
of the signal almost disappear. To measure the transient 
signals, nonlinear optical experiments like four-wave mix- 
ing are usually adopted, in which many-body correlation 
is important sometimes [p6| and may induce new inter- 
esting effects on XS, which, however, is out of the scope 
of our present study. 

Fig. H displays snapshots of the probability distribu- 
tion of the exciton wavepacket resonantly excited by a 
pulse (ilo — E g =- 8.88 meV) with r = 20/uj, corre- 
sponding to the spectral width 0.5 meV. When the NIR- 
excitation pulse is over, the wavepacket performs almost 
perfect periodic swinging along the z— axis with the turn- 
ing points located at about ±zap, except for an overall 
exponential decay due to the ionization. This periodic 
behavior is verified by the time-dependence of the opti- 
cal signal (see the inset). Because of the inter-sideband 
mixing, the wavepacket shape varies drastically within a 



period, in contrast to the unchanged shape of the period- 
ically oscillating wavepacket in the KH potential without 
inter-sideband coupling [Q . All these demonstrate that a 
quasi-bound Floquet state is formed in the strong THz- 
field, and justify our ascribing the suppression of the Fano 
broadening and dephasing rate to the adiabatic XS. 

In summary, the novel excitonic stabilization effect, a 
counterpart of the atomic stabilization effect, has been 
predicted and explored nonperturbatively in a coupled 
system of semiconductors and intense THz-fields by high- 
precise calculation. 
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FIG. 1. Linear absorption spectrum for several THz- field 
strengths indicated by F. The upper-right inset shows, as 
an example, the Fano resonance fitting process for F = 24 
kV/cm, where the squares are the fitted data. The lower-right 
inset shows schematically the sideband structure of the Flo- 
quet-state excitonic spectrum. 

FIG. 2. Fano broadening constant versus the THz-field 
strength. 

FIG. 3. Real-time dependence of the interband polariza- 
tion for various THz-field strength indicated by F, in which 
the curves are offset for clarity. The intensity profile of the 
NIR pulse (the dotted curve) is also shown for comparison. 

FIG. 4. Probability distribution of the NIR-pulse-excited 
wavepacket at several time instant within a period T for 
the THz-field of F = 24 kV/cm. The inset shows the 
semi-logarithmic plots of the intensity of the NIR pulse (dot- 
ted line) and transient interband optical signal (l'P(i)l 2 , solid 
line) as functions of time, and the time snatch (marked by 
the bar) where the snapshots are made. 
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